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A measurement of the decay K'^ — ^ ■K^^'^Vfj,'^ has been performed with the E787 detec- 
tor at Brookhaven National Laboratory. Forty events were observed in the signal region with 
the background expectation of (16.5 ± 2.7) events. The branching ratio was measured to be 
(1.58 ± QAQ{stat.) ± Qm{syst.)) x 10"^ in the kinematic region > 30 MeV and 6*^^ > 20°, 
where E-, is the energy of the emitted photon and 0^^. is the angle between the muon and the 
photon in the K'^ rest frame. The results were consistent with theoretical predictions. 

PACS numbers: 13.20.Eb,12.39.Fe 



I. INTRODUCTION 

A measurement of the decay tt^ v ^i{K ^■i^) 

has been performed using the E787 apparatus [l| at the 
Alternating Gradient Synchrotron (AGS) of Brookhaven 
National Laboratory. This decay is due to radiative ef- 
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fects in the semi-leptonic transition of K'^ Ti^fX^v^ 
decay {K^'^). The decay if^a-y can proceed via inner 
bremsstrahlung (IB) in which a photon is emitted from 
the charged particle in the initial or final state, and 
structure-dependent radiative decay (SD) which involves 
the emission of a photon from the intermediate states in 
the hadronic transition from to tt^. For hadron de- 
cays in the energy region below 1 GeV, the effective-field 
approach based on chiral symmetry, chiral perturbation 
theory (ChPT) % is apphcable Q. At 0{p^) in ChPT, 
a partial branching ratio BR{K^'i^, > 30 MeV and 
Ofj^j > 20°), where Ej is the energy of the emitted photon 
and 9f^^ is the angle between the muon and the photon 
in the rest frame, is predicted to be 2.0 x 10~^ 
The relative size of SD contribution is around 8 %. 

The first search for K^^-y, which was performed in 1973 
Q with decays at rest in a heavy-liquid bubble cham- 
ber, obtained BR{Kf,3^, E-, > 30 MeV) < 6.1 x IQ-^ at 
the 90% confidence level (C.L.). The R- 7r°^-f^7 de- 
cay in flight was studied by the ISTR-I- spectrometer; the 
ratio BR{K^3^)/BR{Kf,3) for the region 30 < E:..^ < 60 
MeV was found to be (4.48 ± 0.68{stat.) ± 0.99{syst.)) x 
lO"'* Based on this value, the Particle Data Group 
cited BR{Kf,3j, 30 < E^ < 60 MeV) as (1.5±0.4) x 10"^ 
% The K+ 7r°/i"'"^'^7 decay at rest was studied 
with the toroidal spectrometer of the E246 and E470 ex- 
periments at the KEK 12-GeV proton synchrotron; the 



BR{Kf,3j, > 30 MeV and > 20°) was found to 
be (2.4 ± 0.5{stat.) ± 0.6{syst.)) x 10"^ [|. 

In this paper, we present results from the E787 ex- 
periment at the AGS. We describe the analysis we 
have performed to observe the ■n^^'^v^^ de- 

cay and the measurement of the partial branching ratio 
BR{Ki,3^, > 30 MeV and 6*^^, > 20°). We compare 
our results to the predictions of ChPT at 0{p'^). The 
BR{K^3j, 30 < Ej < 60 MeV) was also determined and 
can be compared with the previous result. The spectra 
of Kf^s^ observables based on ChPT at 0{p^) 0, IH are 
shown in Fig. [T] and are consistent with the spectra given 
in M. 



II. DETECTOR 

E787 was a rare kaon-decay experiment studying 
K+ -> n+iyp [TTl - [l3j and related decays [H, [lH using 
decays at rest. The data for the K^3^ measurement 
were acquired during the 1998 run of the AGS, using 
kaons of 710 MeV/c incident on the E787 apparatus at 
a rate of about 4 MHz. Two stages of electrostatic par- 
ticle separation in the beam line reduced the pion 
contamination at the E787 apparatus to 25%. The inci- 
dent kaons were detected and identified by a Cerenkov 
counter, multi-wire proportional chambers and a beam 
counter. After being slowed by a BeO degrader, 27% of 
the kaons came to rest in an active stopping target lo- 
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cated at the center of the detector system 
Tesla solenoidal magnet (Fig. [2]) . 

The 12-cm diameter target, which consisted of 0.5-cm 
square plastic scintillating fibers, provided initial track- 
ing of the stopping kaon as well as the charged decay 
products. They passed through a 0.64-cm thick layer of 
plastic scintillation counters surrounding the target (I- 
counters). A 3 x 10~^ radiation length (Xq) thick cylin- 
drical drift-chamber [18] provided tracking information 
for momentum determination. The charged particles lost 
energy in an array of plastic scintillation counters called 
the Range Stack (RS), which provided a measurement of 
kinetic energy and range (equivalent cm of plastic scin- 
tillator) of the charged track. The RS was segmented 
into 24 azimuthal sectors and 21 radial layers, totaling 
one Xq. The RS counters in the first layer (T-counters), 
which were 0.635-cm thick and 52-cm long, defined the 
solid angle acceptance for the charged track in the RS. 
The subsequent RS counters in the second layer (L2- 
counters) and beyond were 1.905-cm thick and 182-cm 
long. The RS counter in the sector and layer where the 
track came to rest was called the "stopping counter" . The 
RS counters were read out by phototubes attached to the 
upstream and downstream ends; the output pulse shapes 
were recorded by 500-MHz sampling transient digitizers 
(TDs) [l^, which provided precise time and energy in- 
formation for the hits in each counter. 

Photons from kaon decays at rest in the target were 
detected by a hermetic calorimeter system surrounding 



the central region. The cylindrical barrel (BL) calorime- 
ter, which was located immediately outside the RS and 
covered two thirds of the solid angle, was used to re- 
construct and measure the three photons from the K^^^y 
decay (i.e. two photons from the tt" and the radiative 
photon). The BL calorimeter consisted of alternating 
layers of lead (0.1-cm thick) and plastic scintillator (0.5- 
cm thick) sheets, and about 29% of the electromagnetic 
shower energy was deposited in the scintillators. The 
BL was segmented azimuthally into 48 sectors and, in 
each sector, four radial groups of 16, 18, 20 and 21 
lead-scintillator layers, respectively with increasing ra- 
dius, formed BL modules totaling 14.3 Xq. The mod- 
ules were 190-cm in length along the beam axis. They 
were read out by phototubes attached to the upstream 
and downstream ends, and the outputs were recorded 
by time-to-digital converters (TDCs) and analogue-to- 
digital converters (ADCs). The two endcap calorime- 
ters 20] , additional "collar" calorimeters for filling minor 
openings along the beam direction, and any active parts 
of the detector not hit by the charged track were used for 
detecting extra particles. 



III. TRIGGER 

The signature of the K^^-y decay was a kaon decay at 
rest with a charged track in the RS and three electromag- 
netic showers in the BL calorimeter. The kinetic energy 
of most of the muons from K^^^ after requiring E-y > 30 
MeV and 61^^ > 20° (Fig. [T]) is less than 85 MeV; thus, 
we used the RS counters in the first six layers for the 
muon measurement. The RS counters in the 7th layer 
and beyond were used to detect extra particles or reject 
photons whose showers started before reaching the BL. 

In the trigger for the Kfj^s^ decay, the "37 trigger", 
a kaon was identified by a coincidence of hits from the 
Cerenkov counter, beam counter and target. To guar- 
antee that the kaon actually decayed at rest, the timing 
of the outgoing particle (via the I-counters) was required 
to be at least 1.5 nsec later than the timing of the in- 
coming kaon (via the Cerenkov counter). A positively 
charged track was required to have a coincidence of the 
hits from the I-counters, a T-counter and an L2-counter 
in the same RS sector (T-2 sector), and from the counters 
in the 3rd or 4th RS layer in the T-2 sector or in either of 
the next two clockwise sectors. The tracks reaching the 
7th or 8th layer were rejected. With such trigger condi- 
tions, only tracks for which the RS stopping counter was 
located in the 3rd to the 6th layer, corresponding to the 
muons with momentum 100-160 MeV/c, were accepted. 

In order to count the number of showers in the BL 
calorimeter online, the analog sum of the outputs from 
8 modules in two adjacent sectors, separately for the 
upstream and downstream ends, were accepted with a 
threshold corresponding to about 5 MeV of visible en- 
ergy per end. The discriminator outputs were OR-ed for 
either end, and then fed into a logic unit. The number 
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FIG. 1: Spectra of -ft"*" — > tv'^ fi'^ u^j observables based on ChPT at 0{p'^): muon momentum (top left), cosine of 9fj,j (top right), 
E-y (bottom left) and neutrino energy (bottom right). The unhatched and hatched histograms represent the distributions before 
and after imposing the conditions E~f > 30 MeV and 9fj,~f > 20°. 




FIG. 2: Schematic side(a) and end(b) views showing the upper half of the E787 detector. C: Cerenkov counter; B4: beam 
counter; I and T; trigger scintillation counters (I-counters and T-counters). Range-stack straw-tube tracking chambers (RSSC) 
were unused in this measurement. 



of showers was required to be equal to three or larger. 
An event was rejected online if the visible energy in the 
endcap calorimeters was more than 20 MeV, the energy 
in the RS sectors outside the region of the charged track 
was more than 10 MeV, or the energy in the RS counters 
from the 9th to 21st layers in any sector was larger than 
the level for minimum ionizing particles. 

The 87 trigger was prescaled by five for taking data 
simultaneously with the trigger for — > n'^h'D. Ap- 
proximately 6 events per slow-extracted beam (2.2-sec 
duration) from the AGS in every 4.2 sec were recorded. 
A total exposure of kaons entering the target available 
for the if measurement was 3.5 x 10^^. A total of 



9.4 X 10^ triggered events were collected. 

Monte Carlo simulation was performed to generate 
samples of K^^sj as well as other kaon decays at rest in 
the E787 detector system. The simulation package de- 
veloped for E787 was comprised of the electromagnetic- 
shower simulation routines from EGS4 [2l[, the routines 
to model other physics processes such as particle de- 
cays and charged-particle interactions in the detector el- 
ements, and the routines to simulate the trigger condi- 
tions. The events in the Monte Carlo samples can be an- 
alyzed with the same analysis codes for real data, except 
that the information on the pulse shapes in the digitiz- 
ers and the information on the outputs from most of the 
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beam instrumentation were not stored. The Monte Carlo 
simulation of E787 is described in [l^ . 



IV. EVENT RECONSTRUCTION 
A. Basic reconstruction 

The charged track was reconstructed and measured 
with a standard algorithm used in the t:^vv anal- 

ysis of E787 described in except that the charged 
particle was assumed to be a muon and all the energy 
deposited in the RS stopping counter was attributed to 
the muon track coming to rest ^ . The initial momentum 
was determined by correcting the momentum measured 
in the drift chamber for the energy loss suffered by the 
muon based on the observed track length in the target. 
The kinetic energy was determined by adding up the en- 
ergy deposits of the muon track in the scintillators of the 
target and the RS, taking account of the energy loss in 
inactive materials such as wrapping and chamber mate- 
rials. The range was calculated from the track length in 
the target and in the RS, and was used for muon iden- 
tification in the if^^a^ analysis. Events were discarded if 
the momentum of the charged particle was larger than 
190 MeV/c. Many of the triggered events were due to 
the — > TT+vr'^ decay for which the 7r+ had interac- 
tions within the RS counters before coming to rest; these 
events were easily rejected because the tt"*" momentum 
(205 MeV/c) was measured with the drift chamber. 

The isolated showers in the BL calorimeter were re- 
constructed and the timing, energy and direction of the 
photons were measured with the same algorithm used 
in the analysis qf^E787 and E949 (the successor experi- 
ment to E787) B2n-2fl for kaon decays with photons in 
the final states [l4l Tla . . The muon track in the RS 
defined the event time reference, and adjacent BL mod- 
ules whose timing was within ±6 usee of the track time 
and whose visible energy was more than 0.2 MeV were 
grouped into a "cluster" . In each module, the energy was 
measured as the geometrical mean of the visible energies 
divided by the calibrated visible fraction (0.292) in the 
upstream and downstream ends. The hit position in the 
end view of the detector was determined by the segmen- 
tation of the modules; the hit position along the beam 
axis {z) was measured from the end-to-end time and en- 
ergy differences. The location of the shower in the BL 
was obtained by taking an energy-weighted average of 
the hit positions in the same cluster and, in conjunction 
with the kaon decay vertex position in the target, the 
azimuthal angle and the polar angle with respect to the 



In the TT"*" track reconstruction, the presence of the additional 4 
MeV deposit in the RS stopping counter due to the 7r+ — >■ 
decay at rest was assumed and was subtracted from the observed 
energy. 



beam axis were determined for each photon. In the K^^^ 
analysis, if two showers were separated by less than 55 
cm and the energy of one of them was less than 25 MeV, 
they were combined to a single shower in order to avoid 
misidentifying a part of an electromagnetic shower to be 
a separate one. Events were discarded unless three or 
four showers were observed in the BL calorimeter. 



B. Primary selection 

Selection criteria ("cuts") were imposed to make sure 
that the trigger conditions were satisfied and the kaon 
in the initial state and the muon and photons in the fi- 
nal state were reconstructed within the fiducial volume 
of the detector system. Furthermore, cuts on the timing 
and energy of the hits recorded in the Cerenkov counter, 
proportional chambers, beam counter and the target, in- 
cluding an offline delayed-coincidence cut which required 
> 2 nsec between the muon time and the kaon time mea- 
sured in the target, were imposed to remove the events 
triggered by kaon decays in flight or by multiple beam 
particles into the detector. From the triggered events, 
1.1 X 10^ events survived and were used in the subse- 
quent analysis. 



V. BACKGROUND SOURCES 

Kaon decays with a single charged track in the region 
100-160 MeV/c and at least one 7r° in the flnal state, i.e., 
/^^3, K+ ^ 7r+^07 {K^2^), K+ ^ n^e+Ve (Kes), K+ ^ 
7r°e+i/e7 (Kea^y) and — > tt+tt^tt*' (-fi'Trs), are potential 
background sources. They are classified as follows: 

"if^a-l": A'ps decays associated with an extra cluster in 
the BL due to accidental hits. 

"i4r^3-2": A'ps decays associated with an extra cluster in 
the BL when the showers due to the two photons 
from TT^ were reconstructed as three showers. 

"^7727" '• -f^Tr27 decay events in which the tt"*" was 
misidentified as a muon or decayed in flight before 
it came to rest in the RS. 

"-^e3(7)" : Kes^y dccays for which the e+ was misidentified 
as a muon, or K^s decays with the e+ misidentifi- 
cation and an extra cluster in the BL or a photon 
due to bremsstrahlung. 

"i^TTrs" : Kj^s decays in which the 7r+ was misidentified as 
a muon or decayed in flight before it came to rest 
in the RS, and at the same time one of the four 
photons from two tt^'s was undetected. 

Background due to multiple beam particles scattering 
into the detector simultaneously was found to be negligi- 
ble. 
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In the subsequent two sections, the cuts to improve 
the signal-to-background and measure the background 
are described. The cuts were developed with the sam- 
ple of real data prescaled by three and the Monte Carlo 
samples. The studies of these samples confirmed that the 
expected number of candidate events in the signal region 
was larger than the background expectation. The cuts 
were then imposed to the entire dataset. 



VI. EVENT SELECTION 
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A. Selection of the muon and three photons 

For muon identification, a cut based on the measured 
range compared to that expected given the measured 
muon momentum was imposed. The measured range 
was also compared to the measured kinetic energy. Fur- 
thermore, to reject e+ tracks in the RS, a set of cuts 
that checked the consistency between the measured en- 
ergy and range in each of the RS counters, the dE/dx 
cuts, were employed. 

For the K^^^ decay, events were discarded unless ex- 
actly three showers with > 5 MeV were observed in the 
BL calorimeter within ±2 nsec of the muon track time 
in the RS. The photon energy should also be sufficiently 
larger than the online energy threshold; since the thresh- 
old had been applied to the analog-sum signals separately 
for the upstream and downstream ends (Section HIT)) , the 
threshold level depended on the z position of the shower 
in the BL. Thus the photon energy was required to be 
greater than 23 (66) MeV for the showers located at the 
edge (center) of the BL. 

In order to remove the events with the photon due to 
bremsstrahlung of the e+ track from the i^es or -ft'e37 
decay within the RS, we calculated the angle between 
the direction of the hit position of the charged track at 
the RS T-countcr from the center of the detector system 
and the direction of each of the photons from the kaon 
decay vertex position. The minimum of the angles should 
be larger than 26° ("e+ bremsstrahlung" cut). 

An extra photon could disappear through inefficiency 
due to very narrow gaps between counters or inactive ma- 
terial. Tight photon-detection requirements, the "photon 
veto" cuts, were imposed on the detector subsystems; hits 
in coincidence with the track time within a few nsec and 
with the energy above a threshold of typically 1 MeV 
were identified as the activity due to an extra photon. 
The photon veto cuts also ensured that only those events 
in which the total photon energy from the K^^-y decay 
was deposited in the BL (and a part of the showers was 
not recorded in any other subsystem) were accepted. 

An extra photon could also go undetected when two 
electromagnetic showers in the BL overlapped each other 
and were reconstructed as a single cluster ("fused clus- 
ter"). In the Monte Carlo studies of K^^^, it was con- 
firmed that the case with the overlap of two photons 
originated from the same ir^ was negligible; thus the 



FIG. 3: Plots of the variables DPSQ vs. CHKSPZ for 
the Kfj,3^ events (left) and the AT^a background events (right) 
generated by the Monte Carlo simulation. 



other case, with the overlap of two photons from different 
7r°'s ("odd combination"), should be considered. We de- 
signed two variables DPSQ and CHKSPZ, described 
below, in which the fused cluster from K^^^ was char- 
acterized. Suppose that one of the three reconstructed 
photons 7i, 72 and 73, e.g. 71, was the fused cluster and 
consisted of two photons 7^ and 7b with the same az- 
imuthal and polar angles and with the energies E^^ — e 
and Ery^ — E^-^ — e, where e was a parameter that varied 
from to E-f-^ . A pair of invariant masses m^^^j and 

and Fi = y/ {111^^72 - ^tto)^ + ("^7b73 - "1^0 )^ 

where m^o is the nominal mass of 7r°, were calculated 
as a function of e. Fi should be small when e agreed 
with the correct value. In the same way, F2 and F3 were 
calculated supposing that 72 and 73 were the fused clus- 
ter, respectively. Finally, the minimum of Fi, F2 and -F3 
with all the possible values of e in each hypothesis, called 
DPSQ in this analysis, was obtained. DPSQ would be 
small if an event was from Kt^^ and one of 71, 72 and 
73 was the fused cluster. On the other hand, the invari- 
ant masses of two out of the three reconstructed photons: 
7717^72, TO7273 and TO7371 were calculated, and the mini- 
mum of |m^j-y2 — TO^o I, ITO7273 — rn^o \ and 1^7371 — rn^a\, 
called CHKSPZ in this analysis, was obtained. If an 
event was from fiT^s with odd combination, the invariant 
mass of other two photons should not be close to 711^^0 
and thereby CHKSPZ would be large. Figure |3] shows 
the DPSQ vs. CHKSPZ plots for the ^^^37 events and 
the if^s background events (after imposing the BL cuts 
and the photon veto cuts) generated by the Monte Carlo 
simulation. To suppress the background, DPSQ was re- 
quired to be more than 30 MeV/c^ and CHKSPZ was 
required to be less than 21 MeV/c^. 

After imposing the cuts in this subsection, 4553 events 
survived. 
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quantity unit resolution 



kinetic energy 


MeV 


0.379 X y/T,, 


momentum 


MeV/c 


\/(0.0227P^)2 - (0.00784P^)4 


azimuthal angle 0^ 


mrad 


17.6 


polar angle 


mrad 


32.2 X cos 6*^ 



TABLE I: Resolutions on the muon observables assumed in 
the kinematic fit to the Kfj,^^ hypothesis. 



B. Kinematic fit 

A least-squares fit was used to improve the resolution 
of the measured it'^a^ quantities. We assumed just one 
undetected particle (neutrino) in the final state. The 
kinematic fit was applied to thirteen observables with 
three constraints. The observables were the kinetic en- 
ergy of the charged particle, the momentum vector of 
the charged particle, and the momentum vectors of three 
photons. The directions of momentum vectors were de- 
fined to be from the kaon decay vertex position. The 
constraints were: 

• The invariant mass of the total energy and momen- 
tum should be equal to the nominal mass of , 

• The kinetic energy and momentum of the charged 
particle should be consistent with a muon hypoth- 
esis, and 

• The invariant mass of one pair of photons should 
be equal to m^o . 

There are three possible pairings of the three photons to 
form the 7r°; the combination which minimized the 
of the global kinematic fit was chosen among the pos- 
sible event topologies. In simulated i^^a-y decays, 86% 
of the events that survived all the selection criteria were 
reconstructed with the correct pairing. 

The kinematic fit was also applied with the assump- 
tion of or i^TrS decay. To identify the K^^ decay, 
the best combination of the two photons from the 7r° was 
chosen as in the case of the fit to the K^^^ assumption, 
and the remaining photon was ignored. Figure Hl^left) 
shows the invariant 77 mass (Mj^) distribution of the 
7r° from the fsT^a events. To identify the if^s decay, the 
events in which four photons were observed in the BL 
calorimeter were used but the photon with the lowest 
energy was ignored; the best combination of two pho- 
tons to form the 7r° was chosen so as to minimize the 
X^- Figure [Upright) shows the AI^^ distribution of the 
7r° from the Ky^s events. The results were used to select 
the events for calibrations and to check the performance 
of Monte Carlo simulation. The resolutions used for the 
muon and the photons in the fit on the K^s^ assump- 
tion, summarized in Table U and Table |lTl were obtained 
from the studies of the if^s events and the Kt^^ events, 
respectively. 

The probability of the kinematic fit on the K^^^ 
assumption, Prob{x\^ ^ ), was required to be more than 
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FIG. 4: M-,^ distributions of the tt" from the K fi,-i{\s,tC) and 
i^TrS (right) events. 



quantity 


unit 


resolution 


energy 


MeV 


1.61 X ^E~, 


azimuthal angle 


j)^ mrad 


39.8 


polar angle Q^, 


mrad 


65.6 X cosS-Y 



TABLE II: Resolutions on the photon observables assumed in 
the kinematic fit to the if^s-y hypothesis. 

10%. In order to suppress the A'^s-l background, the 
events whose x^ probability of the kinematic fit on the 
if^a assumption, Prob{xji^^), was more than 10% were 
discarded. After imposing these probability cuts, 360 
events survived. 



C. Signal region 

The following cuts were applied to the fitted observ- 
ables of the surviving events. The energy of the unde- 
tected particle (neutrino) in the final state, E^, should 
be larger than 60 MeV to suppress the i^,r27 background. 
The invariant mass of the muon and neutrino was calcu- 
lated and, in order to suppress the background events 
with 7r+ decay in flight, the invariant mass was required 
to be larger than 200 MeV/c^. The polar angle of the 
neutrino momentum vector with respect to the beam axis 
should be larger than 26°, to prevent a photon in the fi- 
nal state from escaping along the beam direction. After 
imposing these cuts, 178 events survived. 

Figure [5] shows the cosine of 6*^^ versus plots of 
the events that survived all analysis cuts in real data 
and in Monte Carlo K^^j simulation. The signal region 
was specified with E^ > 30 MeV and 20° < < 60° 
(0.50 < cos(6I^T,) < 0.94). The cut 9^^ < 60° was re- 
quired to improve the proportion of the correct pairing 
to 7r° and to suppress the Kj^^ background. Figure [S] 
and Figure[7]show the Pro6(x|-^^^ ) distributions and the 
M~^j distribution of the 7r°, respectively, from the events 
in the signal region in real data and in Monte Carlo K^s^ 
simulation. Forty events were observed in the signal re- 
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FIG. 5: Cosine of 9^-, versus E~f plots of the events with all analysis cuts imposed in real data (left) and in the sample generated 
by Monte Carlo isT^js^ simulation (right). The box indicates the signal region. 




FIG. 6: Distributions of Prob{x^ 3 ) from the events in the 
signal region in real data (left) and in the sample generated 
by Monte Carlo -ftT^s^ simulation (right). All selection cuts 
except for the cut Pro&(x/f ^,3^ ) > 0.1 in Subsection IVI Bl are 
imposed. 



gion. 



VII. BACKGROUND EXPECTATION 

The K^3-l background was studied with real data from 
the 37 trigger. We used the timing of the photons rel- 
ative to the muon track time, At^. The At^ distribu- 
tions of the two photons to form the 7r° and of the ra- 
diative photon, after imposing all the offline cuts except 
for the timing cut to the photons (|At.^| < 2 nsec in 
Subsection IVIAp . are shown in Figure [51 Figure [HJleft) 
shows the cos{9fj,.y) versus E.y plot of the events with 
a photon in the interval 3 nsec < |Ai^| < 6 nsec in 
real data; the number of events in the signal region in 
Fig. [SJleft) was 11. Assuming that the timing distribu- 
tion of the BL cluster due to accidental hits was con- 



FIG. 7: M^T, distributions of the n from the events in the 
signal region in real data (left) and in the sample generated 
by Monte Carlo -ftT^a^ simulation (right). 



stant, the background level of K^^-l was estimated to be 
ll/(6nsec/4nsec) = 7.3 ± 2.2 events. 

Since it was difficult to isolate samples of the Kf^3-2 and 
_fsr^27 backgrounds from the real data, these backgrounds 
were studied with the samples generated by Monte Carlo 
simulation. The background level of K^j,-2 was estimated 
to be < 0.35 events at the 90% C.L. and was omitted from 
the background estimate. The background level of KT^2'y 
was estimated to be 0.38 ± 0.08 events. 

The Ji'e3(7) background was studied with the real data. 
The dE/dx cuts were inverted to enhance the track 
events, but no event was left after imposing all the other 
analysis cuts. After taking into account the rejection 
of the dE/dx cuts on e"*" tracks as well as the tagging 
efhciency of the inverted dE/dx cuts, the background 
level was estimated to be < 1.1 events at the 90% C.L. 
The background level of i^e3(.^) was omitted from the 
background estimate. 

The _ftr7r3 background was studied with the real data, 
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FIG. 8: At-y distributions of the two photons to form the ir^' 
(left) and of the radiative photon (right) in real data. 



Muon trigger component 

(Section Im)) 
Photon trigger component 

(Section Im)) 
Primary selection 

(Subsection HVBll 
Selection of muon and photons 

(Subsection Ivrai 
Kinematic fit 

(Subsection rvTBll 
Cuts after the kinematic fit 

(Subsection inSl) 
Signal region 

(Subsection inSl) 
Accidental loss 



Samples 


140 


MC 


00609 


MC 


511 


MC, K^2 


220 


MC 


275 


MC 


649 


MC 


351 


MC 


669 


Kf^2, T^scat 



Total acceptance 



4.01 X 10" 



Source 


Background level 


Samples 


a:m3-i 


7.3 ±2.2 


Data 


A-f.3-2 


< 0.35 (90% C.L.) 


MC 




0.38 ±0.08 


MC 




< 1.1 (90% C.L.) 


Data 




8.8 ± 1.6 


Data 


Total 


16.5±2.7 





TABLE III: Expected background levels in the signal region. 
"Data" and "MC" in the rightmost column indicate whether 
real data or Monte Carlo simulation were used for the esti- 
mation, respectively. 



in particular with the muon identification cuts and the 
cuts to detect the extra photon. The studies of photon 
disappearance through the inefficiency of the BL, the in- 
efficiency of the other subsystems, and the fused cluster 
in the BL were made with the BL cuts, the photon veto 
cuts, and the DPSQ and CHKSPZ cuts, respectively. 
Figure right) shows the cos{9^^) versus plot of the 
events after imposing the inverted BL photon-veto cut 
(detection of the activity due to an extra photon in the 
BL calorimeter) in real data. The number of events in 
the signal region in Fig. [HII right) was 45 and, after taking 
into account the rejection of the offiine BL photon-veto 
cut as well as the tagging efficiency of the inverted cut, 
the background level due to the inefficiency of the BL 
was estimated to be 7.5 ± 1.5 events. With the similar 
method, the background levels due to the inefficiency of 
the other subsystems and due to the fused cluster in the 
BL were estimated to be 0.8 ± 0.6 events and 0.5 ± 0.4 
events, respectively. The background level of iC^a was 
8.8 ±1.6 events in total. 

The background levels are summarized in Table Hill In 
total, 16.5 ± 2.7 background events were expected in the 
signal region. 



TABLE IV: Acceptance factors for the A'^s-^ decay in the 
kinematic region E-y > 20 MeV, and the samples used to 
determine them. "MC" in the rightmost column means the 
sample generated by Monte Carlo simulation. "/('^2" and 
"TTscat" mean the data samples of — >■ /i^i'^ decays and 
scattered beam pions, respectively. 



VIII. SENSITIVITY 

The acceptance factors for the selection criteria in this 
measurement (Table IIV|1 were estimated primarily from 
Monte Carlo simulation. We generated the K^^^ sample 
at 0(p4) in ChPT with > 20 MeV, which was set 
to be lower than the offline criteria on the photon en- 
ergy ^. In the simulated events, the kinetic energy and 
momentum of the muon and the energy and z position 
of the photons were smeared with deviates drawn from a 
Gaussian distribution to obtain measured resolutions on 
the observables (Table U and Table |TT| . The acceptance 
of the primary selection to the Cerenkov counter, pro- 
portional chambers, beam counter and the target (Sub- 
section IIVB|1 and the acceptance loss by the accidental 
hits in the detector subsystems, in particular the loss due 
to the veto conditions in the trigger (Section IIIip . were 
measured with data samples of — > pL^v^ decays and 
scattered beam pions, which were simultaneously accu- 
mulated by the calibration triggers. The total acceptance 
was 4.0 X 10~^ and was dominated by the trigger accep- 
tance. 

The single event sensitivity (SES) for if^a-y was de- 
rived from the total acceptance, the total exposure of 



^ Thus, the acceptance factors and the single event sensitivity in 
this section are for the -ff^a-,. decay in the kinematic region E-y > 
20 MeV. In Section llXI we will first obtain the partial branching 
ratio for K^^-^^ in this kinematic region and then convert it into 
BR(Kij,3j, E-y > 30 MeV and e^,-y > 20°) and BR(K^3^, 30 < 
Ej < 60 MeV). 
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FIG. 9: Distributions in cos{9^-f) versus E-, of events selected as the _K'p3-l(left) and iC-n-s (right) background events as described 
in the text. The box indicates the signal region. 



Estimated systematic uncertainties (%) 

BL visible fraction 2.7 

BL z measurement 1.0 

e'^ bremsstrahlung cut 1.0 

Resolution in and 2.8 

Resolution in E.y 1.5 

Resolution in photon z position 1.7 

Monte Carlo statistical uncertainty 1.1 

Total systematic uncertainty 4.9 



TABLE V: Systematic uncertainties in the sensitivity and 
their estimates. 



kaons entering the target (3.5 x 10^^ in Section HIH) . and 
the fraction of kaons entering the target that decayed at 
rest (called the stop efficiency, Fs). Fs was measured 
to be 0.669± 0.014 with the events of the IB component of 
the KT^2'y decay, which were collected with the 87 trigger 
and selected ^ with the method used in the i^7r27 analy- 
sis of E787 [i3, • Thus, the sensitivity for K^s^ was 
normalized to if7r27 * and many systematic uncertainties 
in the acceptance factors canceled. Fg was also measured 



with the Ktts events collected with the 87 trigger, and the 
results were consistent. We obtained SES — 1.07 x 10~^. 
The systematic uncertainties in SES were due to: 

• The BL calibration of the visible fraction and z 
measurement, to which the photon selection (Sub- 
section |vrX| was sensitive, 

• The difference, between real data and Monte Carlo 
simulation, in the distributions of the variables to 
select the muon and three photons, in particular 
the variable on the e"*" bremsstrahlung cut (Sub- 
section |vrX| , and 

• The smearing of the observables in the Monte Carlo 
events when the kinematic fit was applied. 

The systematic uncertainties were studied from the ob- 
served sensitivity variation as each corresponding param- 
eter was varied over the range that had not been ex- 
cluded in the calibration. The statistical uncertainty of 
the Monte Carlo simulation in estimating the acceptance 
was also included. The systematic uncertainties are sum- 
marized in Table fVl By adding them in quadrature, the 
total systematic uncertainty was estimated to be ±4.9%. 



^ The triggered events were further prescaled by eight, and were 
reconstructed with the assumption that the charged track was a 
7r+. After imposing primary cuts, tt"*" identification cuts, photon 
selection cuts and photon veto cuts, the kinematic fit with the 
assumption of the ^^27 decay was appUed. The missing energy 
should be smaller than 141 MeV, and the 7r+ momentum 
should satisfy 140 < < 180 MeV/c. A total of 2425 events 

survived, with the background expectation of 3.6%. The contri- 
bution from the direct photon-emission component of the K.^2-y 
decay was estimated to be 2.2%. 

* The theoretical prediction for the IB component of the -R',r27 
branching ratio in the 7r+ kinetic energy region of 55 to 90 MeV, 
2.61 X 10"'* [2^ . was used in this analysis. The Particle Data 
Group average, (2.75 ± 0.15) X 10"" 0, confirmed this value. 



IX. RESULTS AND CONCLUSIONS 

From the 40 observed events in the signal region, after 
subtracting the expected background (16.5 ± 2.7) we ob- 
tained 28.5 ± 6.9{stat.) K^^^ events with the statistical 
uncertainty of 29%. With SES = (1.07 ± 0.05(s?/si.)) x 
10~^, the partial branching ratio for ^^^37 in the kine- 
matic region E^ > 20 MeV was (2.51 ± 0.74(siat.) ± 
Q.12{syst.)) X 10"^. The kinematic distributions of the 
observed 40 events were compared to the spectra pre- 
dicted from the -ftr^37 decay plus all the background 
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contributions, as shown in Fig. [TOl The x^/degree-of- 
freedom values to evaluate the match between data and 
predictions were 2.7/6, 11.3/8, 12.5/5 and 7.0/9 in the 
distributions of , tt'^ energy, E^, and muon momentum, 
respectively. 

The factors 0.628 and 0.437, estimated from the the- 
oretical K^3^ spectrum at O(p^) in ChPT, are used 
to convert the measured partial branching fraction into 
BR{Ki,3^, E^ > 30 MeV and 6*^^ > 20°) and 
BR{K^3^, 30 < E^ < 60 MeV), respectively. Finally, 
the results of the measurement were: 

BR{Kf,3^, E^ > 30 MeV and 9^,^ > 20°) 
= (1.58 ± 0A6{stat.) ± 0.08{syst.)) x 10"^ 

and 

BRiK^3^, 30<E^ < 60 MeV) 

= (1.10 ± 0.32{stat.) ± 0.05{syst.)) x lO^^. 

The BR{K^,3^, E^ > 30 MeV and > 20°) was con- 
sistent with the theoretical prediction 2. Ox 10~^ 0,1^ as 
well as with the previous measurement (2.4±0.5(stat.) ± 
0.6{syst.))xl0-^ in 3, and the B R{K ,,3^ , 30 < E^ < 60 



MeV) was consistent with the previous measurement 
(1.5±0.4) X 10-5 in 0. The results of the E787 measure- 
ment provide better precision than the previous results. 
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FIG. 10: Spectra of the — s- ty'^^i'^u^j candidates: E~, (top left), tt" energy (top right), Ei, (bottom left) and muon 
momentum (bottom right). In each plot, the circles with error bars represent the distributions of the 40 events observed in 
the signal region. The conditions > 30 MeV and E^ > 60 MeV have been imposed in the analysis. The right-hatched and 
left-hatched histograms represent the expected distributions of the K^s and i^vrs backgrounds, respectively. The unhatched 
histogram represents the distribution of the Monte Carlo Kf^s-y events with the central value of the measured branching ratio 
plus all the background contributions, and should be compared to the circles with error bars. 
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